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Abstract 

The Casimir effect is considered for a wedge with opening angle a, 
with perfectly conducting walls, when the interior region is filled with an 
isotropic and nondispersive medium with permittivity e and permeability 
fj,. The electromagnetic energy-momentum tensor in the bulk is calcu- 
lated, together with the surface stress on the walls. A discussion is given 
on the possibilities for measuring the influence of the medium, via the 
Casimir- Polder force. 
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1 INTRODUCTION 



The wedge geometry is an attractive system to study in connection with the 
Casimir effect, since the geometry is nontrivial enough to exhibit the essentials of 
phenomenological quantum field theory in continuous media, and yet so simple 
that it avoids the formal divergences that so often plague specific calculations 
once curved boundaries are present. An additional bonus from considering a 
system of this kind is that one experiences an interesting formal analogy with 
the theory of a straight cosmic string. 

In this paper we will consider the Casimir theory of a wedge-shaped region of 
opening angle a, when the walls located at angles 6 = and 9 = a are perfectly 
conducting, and the interior region < 6 < a is filled with a homogeneous 
and isotropic dielectric medium of constant permittivity e and permeability fi. 
Figure 1 shows the geometry. The cusp is at the origin. 

The present paper is a generalization of earlier work [1] , in which the interior 
volume was taken to be a vacuum, and is also closely related to Ref. [2]. As 
for references to the earlier literature, we mention those listed in [1] and [2]. 
We may only recall here that our formalism is based upon Schwinger's source 
theory, as developed for the strongly related case of cylindrical geometry in [3] 
and [4], and is related to the general formalism given in Stratton's book [5]. A 
somewhat different approach is followed by Mostepanenko and Trunov in their 
book [6]; their Section 2.3 treats the specific wedge geometry in the vacuum 
case. 

We ought to point out that the formulation of Casimir theory to include 
material properties in the bulk does not imply merely an almost trivial input of 
factors e and fi. The phenomenological electrodynamics is generally quite differ- 
ent from electrodynamics in a vacuum. In particular, the four-momentum of a 
photon in a medium as constructed on basis of Minkowski's energy-momentum 
tensor - or equivalently from the Hamiltonian approach - is spacelike, so it is 
possible to make the electrodynamic field energy negative by means of a Lorentz 
transformation. (Physically, it is precisely properties of this sort that are under- 
lying the recent discussions on the so-called analog models of general relativity; 
cf., for instance, the papers of Leonhardt and Piwnicki [7] and the conference 
report in [8].) 

In the following section we derive the expression for the fundamental dyad F, 
from which the effective products of the fields can be constructed. The two scalar 
Green functions, Fj^ and G^, are determined. In Section 4, the electromagnetic 
energy-momentum tensor is calculated. It is rather remarkable, as shown 
by Eq. (|39|), how the formalism conspires so as to give a very simple result: the 
components of Q^i, reflect the presence of the medium only through a common 
prefactor l/^/e/I- I^i Section 5 we discuss possibilities for measurements, in 
particular, how the deflection of an atomic beam in a medium-filled wedge is 
influenced by the medium. Also, the formal analogy with the theory of cosmic 
strings is briefly commented upon. 

We put h = c = 1, and adopt electromagnetic Heaviside-Lorentz units. 
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2 DYADIC SOLUTION 



Referring to the formalism developed in [1] and [2] for the vacuum fields in the 
bulk, wc can here be brief. We give only the basic definitions, and write down 
formulas when they deviate from the vacuum case. The starting point, as always 
when working with Schwinger's source theory, is the relationship 



E 



(x) = j dx'T{x,x')-'P{x) (1) 



between the electric field E [x) and the polarization source P {x). Here F (x, x') 
is the basic dyad in the formalism; its Fourier transform F (r, r', w) follows from 



/OO 1 
|^e---F(r,r',u;), 



(2) 



with T = t — t'. Maxwell's equations lead to the governing equation 

V X V X r (r, r', w) - efioj'^T (r, r', w) = -iioj^lS (r - r') , (3) 

where 1 is the unit dyad. It is advantageous to introduce a new dyad F' which 
is divergence-free, 

F' = F+il(5(r-r'), V • F' = 0. (4) 

The effective electric and magnetic field products are (in addition to [1-4] cf., 
for instance, also [9]): 

i{E,iv)Ej{v'))^ = F^(r,r',a;), (5) 

i(F,(r)F,(r')L = -^(VxF'xV),, 

= ^-^eifciejmn5fe9^r(„ (r,r',w) . (6) 
fx CO 

Angular brackets mean quantum mechanical expectation values. These effective 
products are to be inserted into the electromagnetic energy-momentum tensor 
{S^v)- The spatial components of are 

Sik = -EiDk - HiBk + ^6ik (E • D + H • B) . (7) 

As F' is divergence-free, it can conveniently be expanded in cylindrical coordi- 
nates. This implies use of the vector spherical harmonics, 

^Im (i^) = [l{l + 1)]"'/' T.Yim (O) (8) 

As mentioned above we let denote the polar angle; the cusp of the wedge 
coincides with the z axis. The boundary conditions are that the electric field is 
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normal, and the magnetic field tangential, at 9 = Q, a (se Fig.l). We let radii r 
and r' correspond to 6, z and d' , z' respectively, and introduce new symbols 
defined by 



m-TT 

Of 



(9) 



with TO a non-negative integer. A lengthy calculation along the lines of [1] and 
[2] leads to the following integral expressions for the spectral dyad F' and its 
double curl: 



m=0 •'-oo 



(V X z) (V' X z) (dm - k^) Fm (r, r') cos cos UmO' 



+— (V X V X 0) (V' X V' X 2-) Gm (r,r') 



(10) 



TO=0 



1 



(V X V X (V' X V' X ?) {dm - k^) Fm {r,r') COSVmOcOSVm^' 



(V X 2) (V X z) {dm - k^) {d'm - k^) Gm {r, r') 

X sinz/„^sinz/TO6''] e*''(''~'''). 



(11) 



Here k G (— oo, oo) is the; axial wave number, the prime on the summation sign 
means that the to = term is taken with half weight, z is the unit vector in the 
z direction, and Fm, Gm are the two scalar Green functions. Further, dm is the 
differential operator 



}_d_ d__K 
dr dr 

The scalar Green functions are explicitly 



Fm (r, r') 
Gm (r,/) 



CO 



ITT 



Sm {r, r') + ^Ju^ {qr<)H,^ {qr>) 



ITT 



Oi^{r,r') + -J.^ {qr<)H,^ {qr>) 



where 



2 2 7 2 

q = e^u) — k , 



(12) 

(13) 
(14) 

(15) 
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and H^^ being ordinary Bessel and Hankel functions of order Um- Further, 
G^*^ the so-called auxiliary Green functions, given by 



-•F,G 



2vn 



, m > 0, 



(16) 



° 2 r> 



(17) 



It ought to be noted that the Green functions are derived on the basis of requir- 
ing boundedness as r — > 0, and outgoing wave conditions as r ^ cx3. 



3 EFFECTIVE FIELD PRODUCTS 

Using Eqs. (^), (^ and ([l^), ( pl| ) we can now calculate the effective field 
product within the wedge, assuming that the two spacetime points x and x' are 
separated. In full generality, all coordinates {i, r, 0, z] would be different from 
{t', r', 0', z'}. However, we assume henceforth t = t' and z = z' , so that x and 
x' are separated only spatially, in the radial and azimuthal directions. 

It is first to be noted that the functions Gq ' do not contribute to the effective 
products. The argument runs similarly to that in [1]. As for the auxiliary Green 
functions we are left only with Gm'~^ with m > 0. 

Now introducing for convenience the operator 



C 



(2^) 



duj 



'• m=0 



we obtain for the diagonal effective products 



i{Eg{r,e)Eg{r',e'))=C 



eg^ rr 



(18) 



(19) 



i{Hg(r,6)Hg{r',6'))=C 



fiq-^ rr 



(20) 



i{Er{r,0)Er{r',e')) 



l^tnfUO^ 1 



jJ^^ (qr^)Hi,^ (qry) ^ sin i^„i6' sin z^m6'' 



(21) 
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i{Hr{r,e)Hr{r',e'))=C 
v'i.eiJ^ 1 



+ 



rr 



-jJ'^m. {Q1'<)Hu^ {qr>) > cos UmO cos UmO' 



(22) 



i{E,{r,e)E4r',e'))=C 



ITT q 

T7 



:{Hjr,e)H4r',0')) 

Y^"^""- '^'^'^<> ^""^ {qr >) cos Vmd cos Urnd' 



(23) 

(24) 
(25) 



These are the diagonal products needed to calculate the normal stresses on the 
surface. There are nondiagonal field products also, (EiEk) and (HiHk), with 
i k. These products do not vanish by themselves, but it is notable that 
the sum e (EiEk) +^ (HiHk), present in the electromagnetic stress tensor, does 
vanish. 

It is also to be noted that the auxiliary Green functions Gm"^ with to > do 
not contribute to the field products. This is because of the relations dmGm'^ = 0, 
as well as 



{dm + eixiJ - e) Fm = u?Gm, 

{dm + enuj^ - k^) Gm = ^Gm- 



(26) 



(27) 



Finally, one may verify by explicit calculation that the Poynting vector vanishes. 
This is as one would expect, under stationary conditions. 

It is to bo emphasized that the differences (r — r') and {6 — 6') are at the 
present stage of the calculation arbitrary; they are not necessarily small. 



4 ENERGY-MOMENTUM TENSOR 

We now let the two points r and r' approach each other, but keep both (r — r') 
and {9 — 9') different from zero. For convenience we write henceforth {E^') 
instead of (Er (r, 9) Er (r', 9'))^^^,, etc. The azimuthal diagonal component of 
{S/j^v) at an arbitrary position r within the wedge is 

{See{r)) = [e (E^^) - e (E^,) + e (E^,) 

+ t,{H',)-n{Hl)+t,{H',)]. (28) 
Inserting the effective products above we get 



{Soe{r)) = \C 



2 I d d d d 

^ rr' 89 d9' dr dr' 



X (gr<) H^^ {qr>)cosum {9 - 9')]^^^, . (29) 
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From this wc have to subtract off the contact term, called (^Sgg (r)). As the 
Casimir effect is caused by the boundaries of the wedge, it follows that the 
contact term has to be evaluated in the absence of any boundaries at all. In 
other words, (5*^^ (r)^ corresponds to a homogeneous dielectric extending over 
all space. Explicit calculation shows that the contact term becomes equal to 
the wedge effective product evaluated at a = it: 



[S'^eir)) = {S^r (r)> 



(30) 



Altough this is a natural result, it could hardly have been written down before- 
hand, without explicit calculation. 

We perform a complex frequency rotation, w iQ, implying 



fc2 



ip. 



(31) 



Then, 



an^ Jo ^0 L I '''' 99 



dr dr' J 



(32) 



where 1^^ and are modified Bessel functions. We go over to polar coordi- 
nates, noting that ^Jepdujdk pdpd(p since p is the radius in the s/eJiQ, k plane. 
Integrating over angles (p from to 7r/2 we obtain 



{See (r)) = 



E 



pdp 



V ^ rn=0 

X I^^ {pr<) (pr>) cos i^m {O - O')] 
We now introduce 

TT 



J_d__d_ _ d_d_ 
rr' 89 89' dr Or' 



P = 



a 



(33) 



(34) 



and assume henceforth that p is an integer. Therewith Vm = mp also becomes 
an integer. This simplifying case is convenient, as it allows us to make use of the 
generalized Graf addition theorem for modified Bessel functions. The important 
formula in our context is given by Eq. (B.6) in [1], and will not be repeated 
here. By means of it, we get 



{See (r)) 



\ / PcJp 

(27r) -^/e/x Jo 



1 8 8 



8 8 



1 p-i 



rr' 89 80' 8r 8r' 



J2 ^0 (z^^") 



m=0 



(35) 
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where 



Rn — 



+ - 2rr' cos ^(61 - 6') + 



1/2 

(36) 



This expression can be processed further, using Eqs. (56) - (60) in [1]. The 
regularized energy-momentum tensor (6^,^ (r)), as defined generally by 

(6^. (r)) = {S,. (r)) - (5°, (r)) , (37) 

then yields for the 99 component 

<«"W> = -™iw(5 + ")(S-0- 

Similar considerations can be carried out for the other components of the energy- 
momentum tensor. If we numerate the components according to (0^,^) = 
{Qrr,Qe9,Qzz, —w) whcrc w is the electromagnetic energy density, we get fi- 
nally 

7r2 



These expressions all vanish for tt = a, as expected. 

The simplicity of the expression ( |39| ) is rather remarkable. The expression 
differs from the corresponding expression in vacuum [1] only through the factor 
y^e/I in the denominator. A value n — > 1 of the refractive index n thus 
causes the energy density in the bulk, as well as the normal stress on the plates, 
to be less than the vacuum value. In the special case of e/x = 1, a case considered 
repeatedly in recent years in various contexts (a so-called " relativistic" medium), 
one ends up with precisely the same energy-momentum tensor as in vacuum. 
There seems to be no simple way to see beforehand why the dependence on the 



properties of the medium should have the special form of Eq. (39) 



5 ON POSSIBILITIES FOR EXPERIMENTS 

We consider first the normal surface density on the lower wall, = 0, as a 
function of the distance r from the cusp. For convenience we will denote this 
force density by a (r) . Since our triplet of basis vectors in the spatial directions 
{r, 9, z} is orthonormal, a (r) must simply be equal to — {Sgg (r)), taken at the 
wall. Thus, in dimensional units. 

The force between the walls is attractive, as expected, and it decreases quickly 
when one moves away from the cusp. The divergence at the cusp is clearly 
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fictitious: the presence of a skin depth 5 in the material forbids us to apply 
the continuum dielectric model at very small distances. Assume copper, for 
instance, for which the conductivity is 6.0 x 10^ {flm)~^. As the most significant 
frequencies contributing to the Casimir force are of the order c/a, where a is 
the local separation between the walls, we obtain J ~ 10 nm if a ~ 1 /zm (the 
typical separation distance) [1]. The expression (|4^) is hardly applicable until 
r becomes as large as about 1 /im. 

Assume for definiteness a very narrow wedge, corresponding to a = lO^^rad 
(0.0057°). A local wall separation of a = 1 /xm then corresponds to r — a/a = 
1 cm. From Eq. ( ^o|) we get 

0.0043 dyn 

a{r = 1 cm) = — — (41) 

y/efi cm"^ 



This is about (Sy/e/l) ^ of the conventional surface force density 0.013 dyn/cm^ 
between two parallel perfectly conducting plates at the same separation as 
above, a = 1 /xm. 

It does not seem to be easy to design a surface force experiment of this kind. 
A more promising possibility might be to consider a variant of the deflection 
experiment of Sukenik et al. [10]. These authors measured the deflection of a 
ground-state atomic beam passing through a wedge-shaped cavity of opening 
angle a ~ 10~'*7'ad, thus of the same order of magnitude as considered above. 
Let us assume the wedge region to be completely filled with a fiuid, having 
material constants e and fi. The transverse force on the beam, causing the 
deflection, will clearly be e and fj, dependent. This transverse force is of course 
the Casimir-Polder rather than the Casimir force in the strict sense, but the 
two kinds of forces are strongly interrelated. 

Let us calculate the interaction energy U (r) for a dipolar particle in the 
wedge as if the particle were at rest, at position r. It is natural to assume that 
retardation effects are not important, so that we can use the material properties 
of the medium at lu = with satisfactory accuracy. Thus, we can use the static 
polarizability a (0) for the particle, and so obtain 

U (r) = -la (0) (E^) = -^a (0) [(Ef.) + {eI) + (E^)] . (42) 

Using Eqs. ([l9|), ( ^ and ( |2^ ) we obtain after some calculation, putting — 0', 

ATTay/eJl e Jq 

CO ' 

2 2 

5 — T^mp {P'r<) Kynp (/or>) cos 2mp9 

- 2Imp {pr< ) Kmp {pr> ) sin^ mpO] . (43) 
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Recall that p = ir/a has been assumed to be an integer. On physical grounds, 
r and r' have to lie close to each other in the expression (p3|). The expression 
may be regularized: if one subtracts off a contact term corresponding to p = 1 
(a = tt), one is left with the wedge-specific contribution, equal to zero in the 
case of a single plane plate in interaction with the dipole. 

The expression ( p3[ ) has the important property that the influence from the 
medium turns up only in the prefactor (y/e/Ie) ^. That is, we can make use 
of the results derived earlier for the case of a vacuum wedge [2]. Defining 
^ = r</r> and taking — 1) to be small, we obtain using Eq. (3.10) in [2]: 



t/(r) 



a(0) 



167r^ 



p4 



2 sin p0 sin p0 



90 



(44) 



Assuming y/ejle > 1 which is usually the case, we thus find a weaker deflection 
of the point dipole toward the nearby wall than in the case of a vacuum cavity. 
In the special case of a plane plate {p — 1), only the first term in (^) survives. 



6 CONCLUSIONS AND FURTHER REMARKS 

Let us first summarize: 

(1) We assumed nondispersive permittivity e and permeability fj, in the 
medium-filled wedge. All calculations were made at zero temperature. The 
electromagnetic energy- momentum tensor (O^i/), after regularization, is given 
by Eq. (|3|). All components of (O^;^) are independent of the polar angle 6 and 
diminish with distance r from the cusp as r~'^. The influence from the medium 
turns up solely in the prefactor .y/ql in the denominator. 



(2) There exists to our knowledge no experiment testing the present kind of 
theory. The most promising possibility seems to be the dielectric variant of the 
experiment of Sukenik et al. [10], measuring the deflection of an atomic beam 
passing through a wedge-shaped cavity of small opening angle a. If the wedge 
is filled with a medium, the transverse force on a dipole is given by the gradient 



of the potential (44), showing that the presence of the medium turns up solely 



in the prefactor y/ejie in the denominator. 
Then a couple of further remarks: 

(3) The formulation of the theory at finite temperatures can be carried out 
in the conventional way, replacing the integral over imaginary frequencies uj by 
a sum over discrete Matsubara frequencies w„ = {innkBT), with n an integer. 
This procedure was shown in [1] in detail for the component {Sgg (r)) , thus 
generalizing the expression (p2|) above to the case of finite temperatures, and will 
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not be further considered here. Again, the important point in our context is the 
apperance of the extra factor y^e/I in the denominator in the energy-momentum 
tensor expression. 

(4) Finally, it is worth noticing that the interesting formal analogy that 
exists between a wedge and a straight cosmic string can be carried over to the 
medium case, only with a slight modification. The line element outside a string 
is 

ds^ = -dt^ + dr"^ + {1 - AGnf r^d0'^ + dz^, (45) 

where G is the gravitational constant and fj, ~ 10^^.g/cm is the string mass per 
unit length (GUT scale). Defining the symbol (3 — {1 — 4G/i)" , we can write 
the electromagnetic energy- momentum tensor as [11] 

(e^.) = ^20^ + 11) - 1) diag (1, -3, 1, 1) . (46) 

Comparison with ( |39| ) shows that, apart from the extra prefactor y^e/I in the 
denominator of that equation, there is complete analogy between the two cases 
if the gravitational quantity f3 is identified with the wedge quantity n/a. 
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